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I.  INTRODUCTION 


This  report  deals  with  an  analysis  of  the  high  frequency  (h.f.) 
radiation  from  slots  in  planar  and  cylindrical,  perfectly-conducting 
surfaces  which  are  partly  covered  with  an  impedance  surface  patch.  The 
present  study  is  applicable  to  the  problem  of  controlling  the  polari- 
zation and  pattern  shape  of  the  electromagnetic  radiation  from  a flush 
mounted  airborne  antenna  for  satellite  communication  purposes.  In  this 
application,  a slot  antenna  is  flush  mounted  in  the  aircraft  fuselage, 
and  one  is  interested  in  the  radiation  pattern  in  the  roll  plane  con- 
taining the  antenna.  Thus,  the  problem  is  essentially  two-dimensional 
(2-D)  in  nature,  and  the  aircraft  fuselage  may  be  approximately  modeled 
by  a circular  cylinder  in  the  roll  plane.  The  aircraft  wing  structure 
is  ignored  in  this  initial  study;  however,  it's  effects  can  be  readily 
incorporated  in  the  future  via  the  geometrical  theory  of  diffraction 
(GTD)  as  will  be  discussed  later  on.  In  particular,  the  behavior  of 
the  radiated  fields  in  the  vicinity  of  the  horizon  (or  the  shadow  boundary) 
associated  with  the  antenna  are  of  prime  interest  in  the  satellite  communi- 
cation type  application.  Thus,  it  is  desirable  for  the  airborne  antenna 
to  have  a gain  approximately  that  of  the  pattern  maximum  near  the  horizon; 
such  a desired  high  gain  may  be  achieved  by  an  impedance  surface  type  of 
loading  around  the  slot. 

Figures  1 and  2 illustrate  the  antenna  configurations  of  interest 
which  are  analyzed  in  this  report.  The  slot  is  present  in  the  perfectly- 
conducting  surface  and  the  perfectly-conducting  surface  is  coated  with 
a sufficiently  long  impedance  surface  patch  which  also  covers  the  slot 
as  indicated  in  Figures  1 and  2.  With  a proper  choice  of  the  impedance, 
the  slot  may  be  allowed  to  radiate  like  an  end-fire  surface  wave  type 
antenna  near  the  horizon;  it  is  this  end-fire  effect  which  may  be 
employed  to  increase  the  gain  near  the  horizon.  By  a proper  choice  of 
the  impedance,  it  is  implied  that  the  type  of  impedance  chosen  is  such 
that  it  supports  a bound  surface  wave  mode  on  the  planar  impedance 
surface  of  Figure  1;  whereas,  it  supports  an  Elliot  type  mode  [1] 
on  the  curved  impedance  surface  of  Figure  2.  The  Elliot  mode  propagates 
along  the  curved  impedance  surface  with  negligible  attenuation  or  leakage, 
and  it  corresponds  in  the  limiting  case  of  the  infinite  radius  of 
curvature  to  the  bound  surface  wave  mode  which  exists  on  the  planar  im- 
pedance surface.  It  is  noted  that  these  surface  wave  type  modes  are 
excited  by  an  axial  slot  in  the  configurations  of  Figures  1 and  2,  if 
the  surface  impedance  is  inductive.  On  the  other  hand,  the  circum- 
ferential slot  will  excite  these  modes  only  if  the  surface  impedance  is 
capacitive.  These  surface  wave  modes  on  the  impedance  surface  diffract 
from  the  ends  of  the  impedance  patch  (as  a result  of  the  discontinuity 
in  surface  impedance  there),  thereby  producing  an  end-fire  effect. 

While  the  cylindrical  antenna  configuration  of  Figure  2 is  important 
in  this  study  because  it  models  an  aircraft  fuselage  in  the  roll  plane, 
the  planar  configuration  of  Figure  1 is  also  analyzed  in  this  report 
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Figure  1.  A thin,  2-0  slot  in  an  infinitely  long,  planar, 
perfectly-conducting  structure  with  a surface 
impedance  patch  of  length  L.  The  value  of  the 
surface  impedance  is  Z5. 


Figure  2.  A thin,  2-0  slot  in  a perfectly-conducting  circular 
cylinder  with  a surface  impedance  patch  of  arc 
length,  L.  The  value  of  the  surface  impedance  is  Zj 


f 
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because  it  is  a more  basic  geometrical  configuration  involving  a perfectly-  j 

conducting  structure  with  a surface  impedance  patch,  and  because  a com-  | 

parison  of  the  patterns  for  the  configurations  in  Figures  1 and  2 allows  ' 

one  to  ascertain  the  effects  of  surface  curvatures  on  the  patterns 
especially  in  the  region  of  space  in  which  the  slot  is  directly  visible. 

In  the  present  analysis,  the  impedance  surface  patch  is  used  to 
approximately  simulate  the  effects  of  a thin,  uniform  dielectric  or  ferrite 
cover  or  properly  designed  corrugated  surface  of  finite  extent  on  the 
perfectly-conducting  planar  or  cylindrical  antenna  structure.  The  slot 
in  the  perfectly-conducting  structure  is  covered  by  the  dielectric  or 
ferrite  material,  and  it  is  assumed  that  the  electric  field  in  the  slot 
aperture  is  known,  so  that  one  also  knows  the  equivalent  magnetic  current 
in  the  aperture.  It  is  noted  that  the  equivalent  magnetic  current  in 
the  aperture  of  the  slot  is  z-directed  for  the  2-D  axial  slot,  whereas 
it  is  T-directed  for  the  2-D  circumferential  slot.  Hence,  the  axial 
slot  may  be  referred  to  as  a z-directed  slot,  and  the  circumferential  slot 
may  be  likewise  referred  to  as  a x-directed  slot.  The  radiation  from  a 
slot  of  finite  width  may  be  obtained  by  quantizing  the  source  distribution 
in  the  aperture,  and  by  superposing  the  fields  radiated  by  each  of  the 
quantized  sources.  An  equivalent  quantized  source  in  this  case  is  a 
z-directed  magnetic  line  source  for  the  axial  slot,  whereas,  it  is  a 
T-directed  magnetic  line  dipole  for  the  circumferential  slot.  These 
equivalent  quantized  magnetic  currents  radiate  in  the  presence  of  the 
perfectly-conducting  surface  (which  covers  the  aperture  as  well),  and  in 
the  presence  of  the  dielectric  or  ferrite  material  which  covers  the 
perfectly-conducting  surface.  The  strengths  of  the  equivalent  sources  are 
weighted  according  to  the  source  distribution  in  the  aperture.  In  the 
present  report,  only  the  radiation  from  a magnetic  line  source  and  a 
magnetic  line  dipole  will  be  analyzed  since  the  radiation  from  a slot  of 
finite  width  may  be  readily  obtained  from  this  analysis  by  superposition. 

In  the  impedance  surface  approximation  for  the  truncated  dielectric  or 
ferrite  cover  on  the  perfectly-conducting  surface,  the  equivalent  problem 
to  be  analyzed  consists  of  an  equivalent  magnetic  line  or  line  dipole 
source  on  the  impedance  surface  patch  which  covers  the  perfectly-conducting 
planar  or  cylindrical  surface. 

In  the  absence  of  any  surface  impedance  patch,  the  axial  slot  in 
an  electrically  large  perfectly-conducting  convex  cylinder  would  radiate 
fields  which  are  substantially  stronger  in  the  vicinity  of  the  horizon 
than  those  radiated  by  a circumferentially  directed  slot  at  the  same 
location.  In  fact,  the  difference  in  the  power  radiated  by  these  two 
slots  could  be  as  much  as  20  dB  or  so  in  the  neighborhood  of  the  horizon; 
and,  in  the  deep  shadow  region  this  difference  continues  to  increase 
significantly.  The  axially  directed  slot  in  a convex  cylinder  radiates 
a (ji-directed  electric^field;  on  the  other  hand,  the  circumferential 
slot  radiates  only  a z-directed  electric  field  in  the  roll  plane.  Thus, 
the  electric  field  components  radiated  by  the  two  slots  are  polarized 
orthogonal  to  each  other.  It  is  observed  that  an  alternate  geometry  that 
could  prove  interesting  is  the  electric  line  source  placed  on  an  impedance 
surface.  This  source  would  in  general  excite  a loosely  bound  surface  wave 
more  strongly  than  a magnetic  line  dipole  source. 
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In  some  satellite  communication  type  airborne  antenna  applications, 
it  is  commonly  desirable  for  the  axial  and  circumferential  slots  to  be 
used  in  combination  as  a crossed  slot  antenna  for  the  purposes  of 
transmitting  or  receiving  a circularly  polarized  (CP)  wave  near  the 
horizon.  However,  the  20  dB  or  so  difference  in  the  power  radiated 
individually  by  the  two  slots  makes  the  crossed-slot  configuration  not 
too  well  receptive  to  a CP  wave  (or  by  reciprocity,  the  crossed-slot 
does  not  effectively  radiate  a CP  wave)  near  the  horizon.  The  present 
analysis  could  be  used  to  ascertain  if  one  could  enhance  the  circumfer- 
ential slot  radiation  near  the  horizon  via  the  end-fire  effect  by  employing 
a capacitive  type  surface  impedance  patch  around  the  slot  (as  in  Figure  2), 
while  simultaneously  not  degrading  the  level  of  the  radiation  pattern  of 
the  axial  slot  component  in  the  crossed  slot  configuration. 

The  antenna  problems  in  Figures  1 and  2 are  analyzed  in  this  report 
via  a combination  of  the  method  of  moments  (MM)  and  the  GTD;  this 
combination  will  be  referred  to  as  the  MM-GTD  method.  An  integral 
equation  is  formulated  for  the  "scattering  problem"  in  which  the  unknown 
tangential  magnetic  field  induced  on  the  surface  of  Figures  1 and  2 by 
a distant  source  needs  to  be  evaluated  only  over  the  surface  impedance 
patch.  In  this  scattering  problem,  the  equivalent  magnetic  line  or  line 
dipole  source  (corresponding  to  each  of  the  quantized  equivalent  currents 
in  the  slot  antenna  problem)  is  turned  oft. 

Since  the  unknown  needs  to  be  evaluated  only  over  the  surface 
impedance  patch,  . le  is  able  to  treat  electrically  large  cylinders  and 
infinite  planar  structures  (of  Figures  2 and  1)  far  wore  efficiently 
than  by  MM  solutions  of  conventional  integral  equations  employing  the 
free  space  Green's  function,  in  which  the  unknown  must  be  evaluated  over 
the  perfectly  conducting  portion  of  the  surface  as  well.  The  source  term 
in  this  integral  equation  represents  the  tangential  magnetic  field  induced 
(by  the  distant  source)  on  the  same  surface,  but  in  the  absence  of  the 
surface  impedance  loading,  i.e.,  on  the  "unperturbed"  or  completely 
perfectly  conducting  surface.  The  kernel  of  this  integral  equation  repre- 
sents the  tangential  magnetic  field  which  is  excited  on  the  "unperturbed" 
surface  due  to  an  appropriate  magnetic  line  source  on  the  same  surface. 
Since  the  source  term  and  the  kernel  of  this  integral  equation  correspond 
to  the  unperturbed  or  the  perfectly-conducting  surface  without  the  surface 
impedance  patch,  they  are  readily  calculated  via  the  special  Green's 
functions  for  the  planar  and  circular  perfectly  conducting  surfaces. 
Furthermore,  when  the  asymptotic  high  frequency  (or  the  GTD)  approximations 
are  employed  for  the  cylinder  Green's  functions,  one  is  able  to  analyze 
not  only  electrically  large  circular  cylinders  more  efficiently,  but  also 
arbitrary,  convex,  perfectly-conducting  cylinders  with  an  impedance  patch. 

The  unknown  surface  field  is  solved  for  via  the  usual  MM  procedure 
over  just  the  surface  impedance  patch  region.  Once  this  unknown  induced 
surface  magnetic  field  is  determined  in  the  "scattering"  problems,  it 
also  directly  furnishes  the  radiation  patterns  for  the  "antenna"  con- 
figurations of  Figures  1 and  2 via  reciprocity. 
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The  analytical  details  are  presented  in  Section  II  where  it  is 
indicated  that  the  circumferential  slot  case  presents  certain  numerical 
difficulties;  a few  approaches  for  circumventing  these  difficulties  are 
proposed,  and  some  preliminary  numerical  results  for  both  the  axial  and 
the  circumferential  slot  cases  are  oresented  in  Section  III. 

II.  ANALYTICAL  FORMULATION 

The  solutions  to  the  antenna  problems  of  Figures  1 and  2 are 
formulated  indirectly  in  terms  of  the  solutions  to  their  corresponding 
scattering  problems  via  reciprocity,  as  indicated  below.  This  formu- 
lation is  presented  for  the  3-D  case,  and  it  is  subsequently  specialized 
to  the  2-D  case  which  is  of  interest. 

Consider  a_surface  S on  which  is  placed  an_infini tesimal  magnetic 
current  moment  Mg  at  Q'  as  shown  in  Figure  3.  is  tangential  to 


Figure  3.  An  illustration  of  the  antenna  and 
corresponci ng  scattering  problems. 


i at 


Q' 


_ _ The  equivalent  source  Ma  generates  the  electromagnetic  fields 

Ea,  ha  in  the  presence  of  S,  and  it  constitutes  the  antenna  problem 
similar  to  that_in  Figures  1 and  2.  Consider  another  infinitesimal 
current  source  F at  point  P wnich  is  external  to  S;  the  ve£tcr  - could 
represent  either  an  infinitesimal  electric  current  moment  J,  or  an 
infinitesimal  magnetic  current  moment  M.  The  source  Ma  in  the  antenna 
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problem  generates  Eg,  Hg  with  the  source  F jturned  off.  On  the  other  hand, 
F generates  the  electromagnetic  fields  Ft,,  Hb  in  the  presence  of  S,  but 
with  the  source  Mg  turned  off;  thus  F constitutes  the  vector  source  for 
the  corresponding  scattering  problem.  All  of  the  above  fields  and  sources 
satisfy  an  eJ“^  time  dependence  which  is  assumed  and  suppressed  in  the 
present  analysis. 

Let  the  surface  S and  the  surface  £ at  infinity  bound  the  volume  V. 
The  electromagnetic  fields  satisfy  Maxwell's  equations  in  V,  and  the 
radiation  condition  on  I.  In  the  present  analysis,  the  electromagnetic 
fields  are  chosen  to  satisfy  th^  im£edance  boundary  condition  on  S.  To 
distinguish  between  the  fields  Eb,  Hb  due  to  the  source  F which  is 
either  J or  M,_one  lets  E^,  H|  be  the  fields  due  to  J,  and  E^,  be  the 
fields  due  to  M,  respectively.  The  Maxwell's  equations  are 


7xE  = 

d 

- - M^; 

= ioeE 
a a 

(la;lb) 

VxE^  = 

- ; 

7xH^  = >eE^+J 

, if  F=J 

(2a;2b) 

and 


VxEj  = 


- jwuHj  - M ; 


VxTiJ  = joeE^  , if  F=M 


(3a;3b) 


The  impedance  boundary  condition  on  S dictates  that 
F,  - (F, ‘n)?!  = nxF  , on  S; 

a a S a 


(4a) 


and  likewise 


re,ni  /Fe,ni 


- (Ef’‘"-n)n  = nxH 


, on  S, 


(4b) 


where  Zg  is  a complex  valued  scaler  function  which  denotes  the  value  of 
surface  impedance;  the  value  of  Zj  could  vanish  over  a part  of  S 
(Zg  = 0 implies  a perfectly  conducting  surface).  From  Equations  (1)  and 
(2) , one  may  obtain 


’ - K-^  - Ea-7 


if  F=J. 


(5) 


Similarly,  one  obtains  from  Equations  (1)  and  (3)  the  following 


J 


Integration  of  both  sides  of  Equations  (5)  and  (6)  over  the  volume  V,  and 
an  application  of  the  divergence  theorem  yields: 

S V 

and 


-ff  - ^*“al  • Sds  » - - Hj-fT]dv  ; if  T-IT.  (8) 

s V 


The  contributions  to  the  integrals  from  the  surface  l at  infinity  vanish 
via  the  radiation  condition.  The  surface  integrals  in  Equations  (7)  and 
(8)  vanish  via  the  impedance  boundary  condition  in  Equation  (4)  as 
follows : 


(Z^nxnxH^  )-H^  • (Z^nxnxH^  =0 . 
The  infinitesimal  sources  M and  T may  be  represented  as 


(9) 


^=V5([R-Rq,1) 


F = 


J=/U(S(|R-Rpi)  ; if  F=J 


M=mv5(lR-Rp!)  ; if  r=M 


_ . (10a;10b) 


where  the  amplitudes  m^,  / and  m are  assumed  known,  an_d  6 denotes  the 
Dirac  delta  function  in  terms  of  the  position  vector  R=xx+yy+zz  in  the 
3-D  coordinate  space.  Thus,  Equations  (7)  and  (8)  become  via  Equations 
(9)  and  (10),  the  following: 


hJ(Q' iP)-t  m^  = - E^(p1Q')-G  j ; if  r=T 

and 

^(Q'lP)-tm^=  H3(PiQ')-vm  ; if  r=M 


(11) 


(12) 


The  notation,  Hb’^(Q'|P)  implies  that  it  is  the  field  produced  by 
_the  source  at  P,  and  evaluated  or  sampled  at  Q'.  Similarly,  Ea(p';Q'), 
8a(P|Q')  are  the  fields  due  to  the  source  at  Q',  and  they  are  evaluated 
or  sampled  at  P.  The  quantities  on  the  LHS  of  Equations  (11)  and  (12) 
may  be  interpreted  as 
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where  the  point  Q"  is  external  to  S.  For  the  sake  of  convenience  one 
lets  y = l_(if  F=0}_,^and  m = _]_  (if  r=Fr) ; it  is  then  clear  that  by 
orienting  Jf_uJ  or  M=vM  along  E,  or  Hg,  respectively  one  can  calculate 
the  fields  Eg,  or  of  the  radiation  problem  which  are  originally 
of  interest,  if  the  induced,  tangential  magnetic  fields  H^(Q'lP)*t»  or 
H]5(Q'lP)’t  in  the  scattering  problem  are  known. 

One  may  directly  employ  the  three-dimensional  (3-0)  results 
developed  above  to  the  2-D  configurations  of  interest  in  which  the 
geometry,  the  source,  and  the  fields  have  no  variation  along  the 
^-direction.  In  the  2-D  case,  the  equivalent  magnetic  current  moment 
Ma  becomes  a magnetic  line  source.  For  the  infinitesimally  thin  axial 
(or  z-directed)  slot  case,  the  equivalent  magnetic  line  source  possesses 
magnetic  current  elements  which  are  z-directed,  whereas  in  the  infinitesimal 
circumferential  (or  t-directed)  slot  case,  the  equivalent  magnetic  line 
source  possesses  magnetic  current  elements  which  are  T-directed  so  that  the 
equivalent  source  is  referred  to  as  a magnetic  line  dipole  in  this  case. 
Thus,  in  th_e_2-D  infinitesimal  axial  slot  case  ^ is  given  by 
Ma=z  ma5( 1p-pqi 1 ) , with  t=z,  and  one  chooses  F=M  with  v=z  since  the  axial 
slot  generates  only  a z directed  ITg  in  the  2-D  problem.  With  m=l.  Equation 
(12)  reduces  in  the  2-D  case  to 

H',(p1Q')’Z  = Hk(Q'  1p)‘Z  ni  ; for  the  axial  slot  case,  (13) 

where  it  is  noted  tha;t  jFor  F=M=z6  ( iF-Fp  1 ) , lib  is  also  z-directed  in  the 
2-D  problem.  The  5 ( jp-pgiQ^pj_)  implies  a 2-0  delta  function  in  the  2-D, 
(x,y)  coordinate  space  with  p=xx+yy  in  contrast  to  the  3-D  delta 
functions  in  the  (x,y,z)  coordinate  space  as  in  Equations  (10a;10b). 

In  the  2-D  infini_tesimal  ci_r£umferential  slot  case,  is  a magnetic 
line  dipole  given  by  Ma=T  ma5(  Ip'-pqi  1 ) with  t=T  being  the  unit  tangent 
vector  to  _yie  2-£  boundary  S^in  the  (x,y)  plane  as  in  Figure  3,  and_one 
chooses  F=J=z6(  [p-ppl ) with  u=z  and  y=l , since  now  generates  an  E^ 
which  has  only  a z-component  in  the  2-D  problem.  Also,  r|  generated  by 
this  T is  entirely  z-directed  in  the  2-D  case.  Consequently,  Equation  (11) 
becomes 

r (p|Q')*z  = - Hf(Q'|P)-T  m,  ; for  the  circumferential 
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As  a next  step,  the  solution  to  the  quantities^Hb(Q' ;P) ‘Z  of 
Equation  (13)  for  the  2-0  axial  slot,  and  H^(Q'iP)*t  of  Equation  (14) 
for  the  2-0  circumferential  slot,  respectively,  will  be  formulated  in 
terms  of  integral  equations  for  the  scattering  problem;  and  these 
integral  equations  will  be  numerically  evaluated  via  tjie  method  of 
moments  (MM),  thereby  also  providing  the  solution  for  Eg,  Ha  of  the 
antenna  problems  via  Equation  (13)  and  Equation  (14)  for  the  two  dif- 
ferent types  of  slots.  ^The  scattering  problem  corresponding  to  the 
evaluation  of  H^(Q‘lP)'Z  on  S will  be  referred  to  as  the  TZz  case; 
whereas,  that  corresponding  to  the  evaluation  of  H|(Q‘iP)'T  on  S will  be 
referred  to  as  the  TMz  case.  The  formulation  leading  to  the  integral 
equations  is  first  given  for  the  3-0  case,  and  it  is  later  specialized 
to  the  2-0  scattering  problems  of  interest. 


(a)  ACTUAL  OR  PERTURBED 
SCATTERING  PROBLEM 


(b)  UNPERTURBED  PROBLEM 


Figure  4.  The  actual  scattering  problem,  and 
a related  unperturbed  problem. 

Figure  4(a)  illustrates  the  same  scattering  problem  as  in  Figure  3, 
except  that  the  surface  impedance  Zs  is  now  specifically  taken  to  be 
nonzero  on  the  portion  S],  whereas,  it  is  zero  everywher9_cn  the 
_r9main_ing  portion  S2  of  the  surface  3.  Again  the  source  F may  be  either 
^ or  M,  and  J radiates  the  fields  ?§,  ri|  in  the  pt-esence  of  S,  likewise, 
M radiates  1^,  as  before. 
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The  scattering  problem  in  Figure  4(a)  will  be  referred  to  as  the 
perturbed  problem  due  to  the  fact  that  the  presence  of  a nonzero  2$ 
over  S]  may  be  viewed  as  a perturbation  of  the  perfectly  conducting 
surface  which  would  otherwise  have  existed  over  the  entire  surface  S. 

Figure  4(b)  illustrates  the  unperturbed  problem  in  which  the 
entire  surface  S is  perfectly  conducting.  The  excitation  chosen  for  the 
unperturbed  problem  in  Figure  4(b)  is  an  infinitesimal  magnetic  current 
moment  TTj  external  to  S whi^h  generates  the  fields  Eqi  Hq  in  the 
presence  of  S.  The  fields  Eq.  Hq  satisfy  Maxwell's  equations 

VxFq  = - jww  Hq-Mj  ; VxHq  = juE  Tq  . (15a; 15b) 

Furthermore,  Tq,  Fo—SsTisfy  the  radiation  condition  over  the  surface 
I at  infinity,  and  Eq  satisfies  the  following  boundary  condition  on  S: 

n X = 0,  on  S.  (16) 

As  before  Equations  (2)  and  (15)  yield 

Similarly  Equations  (3)  and  (15)  yield 

V-Cr^xfiJ  - eJxH-q]  = - H^-M^  + Hq-M  ; if  F=rT  . (18) 

An  application  of  the  divergence  theorem  to  the  LriS  of  Equations 
(17)  and  (18),  and  the  use  of  the  radiation  condition  on  Z,  the  im- 
pedance boundary  conditions  (Equations  (4a;4b))  on  Si  with  Zs=0  over  S2, 
and  the  use  of  the  boundary  condition  in  Equation  (16)  yields: 

jJn’^-EZgnxnxH^lds  = - H^) 

S^  V V 

and 

-[Z^nxnxHjlds  = - JJJ' Hj-M^dv  + JJ'f  -Mdv  ; if  . (20) 

S^  V V 
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F is  as_  defined  previously  in  Equation  (10b),  and  y=l , m=l  for  convenience. 
Also,  Rj  is  chosen  to  be  a unit  amplitude  test  source. 
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= t 6(lR-Rp. I).  (21) 

Thus,  Equations  (19)  and  (20)  respectively  become 

(^lP')*[Zs(Q)BxnxH^(Q|P)]ds  = - H^(P'iP)*T-r  (PlP')-G  ; 

^1 

‘ if  F=J  , (22) 

and 

j[jrH^(QlP')-[Z^(Q)nxnxHj(QlP)]ds  = - Hj(P‘ 1 P)  •t+TT^(P  1 P' ) • v ; 

^1 

' if  F=M.  (23) 

In  Equations  (22)  and  (23),  represents  a point  on  S-j  and  it  is  there- 
fore a variable  of  integration;  also  the  notation  for  the  field  quantities 
in  terms  of  the  points  Q,  P,  and  P'  has  the  same  meaning  as  before  (see 
Equations  (11)  and  (12)). 

One  may  now  specialize  Equations  (22)  and  (23)  to  the  2-0  TE^,  and 
TMz  scattering  problems  of  inter,e^t  by  assuming  no  variation  in  tne 
^-direction.  In  the  TEz  case,  F=M  with  v=z;_and  in  the  TMz  case, 

F=T  with  u=z.  F^rtherm£re_,  the  test  source  My  in  the  unperturbed  problem 
is  chosen  to  be  My=z6( [p-ppi 1 ) with  t=Z  in  the  TEz  case,  thereby  consti- 
tuting a unit  amplitude  magnetic  line  source  at  the  other  hand, 

My^is  chosen  to  be  a magnetic  line  dipole  %=t5(  Ip-ppi  |)  at  P'  such  that 
T=t  (a  unit  tangent  vector  to  the  2-D  boundary  S illustrated  in  Figure  4) 
as  the  point  P'  moves  to  the  point  Q'  on  the  boundary  S.  In  fact,  the 
required  integral  equations  result  when  P'  moves  to  Q'  on  S in 
Equations  (22)  and  (23);  thus, 

/HQ^(2!iQ')-[Z^(Q)nxnxRj(Q|P)]dA  = - \P) \Q' ) -1  ; 

^1 

for  F=M,  or  TE^  case,  (24) 
and 


11 


^(QlQ-)-[Z^(Q)fixnxH^(Q)lP)3di  = - H^(Q' 1 P) iQ' ) 

Si 

for  r=J,  or  TM^  case.  (25) 

In  order  to  distinguish  between  the  fields  due  to 

Mj=z5 ( 1 ) for  the  TE^  case,  or  Mj=t6(  U-Jg, | ) for  the  TM^  case, 

one  defines  as  being  produced  by  Mj=26(1F-Fqi 1)  for  the  TE^ 

case;  and  likewise,  one  defines  as  being  produced  by 

R^=to( jF-Fqi I ) for  the  TM^  case,  as  indicated  in  Equations  (24)  and  (25). 

It  is  noted  that  ds=dJid2  where  di  is  an  incremental  arc  length  along  the 
boundary  curve  S'] ‘for  the  2-D  case,  and  the  integral  over  dz  is  sup- 
pressed in  the  2-D  case.  In  a given  coordinate  system,  the  direction 
of  integration  along  dJi  is  chosen  to  increase  in  the  positive  sense  such 
that  ds=didz>0  with  dz  defined  as  being  in  the  increasing  z-direction. 

To  be  specific,  let  dl^xldTl  so  that  t may  now  represent  the 
increasing  coordinate  direction  for  di.  The  latter  aspect  will  be 
clarified  further  in  the  specific  applications  to_planar  and^cylindrical 
configurations  of  Figures  1 and  2.  Since  fTj  and  M are  both  z-directed 
in  the  TE2  case,  the  fields  Hqz  and  HE_are  also  entirely  z-dj_rected 
in  Equation  (24).  On  the  other  hanci,  Mj  is  t-directed,  and  J is 
z-directed  in  the  TMz  case;  thus,  Eqt  is  entirely  z-directed.  Further- 
more, itjnay  be_shown  by  applying  the  reciprocity  theorem  to  the  pair  of 
sources  Mj  and  F in  which  both  Mj  and  F now  radiate  in  the  presence  of 
the  unperturbed,  or  the  perfectly  conducting  surface  S (corresponding 
jto  Figure  4(b))  that  the  field  Fioz(PlQ')  in  Equation  (24),  and  the  field 
Eox(PlQ')  in  Equation  (25)  are  related  to  the  fields  FJ[j(Q'iP)  and 

H^^(Q' [P)  as  follows : 


Hoz(PlQ')-z  = ’ (26) 

and 

Eq^(P|Q')-z  = - H^^(Q' IP)-^  . (27) 


The  field  H^b(Q' 1 P)  is  generated  by  F=M=z6( Ip-ppj ) at  P in  the 
presence  of  the^unperturbed  surface  S (Figure  4(b))  for  the  TEz  case,  and_ 
it  is  entirely  z-directed;  whereas,  H^5(Q'|P)  is  generated  by  F=T=z6(ip-pp 
at  P in  the  presence  of  the  unperturbed  surface  S (Figure  4(b))  for 
the  TMz  case.  The  results  in  Equations  (26)  and  (27)  are  similar  to 
those  in  Equations  (13)  and  (14),  respectively;  in  fact,  they  may  be 
derived  in  a manner  identical  to  that  employed  to  arrive  at  Equations 
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01)  and  (12)  (or  (13)  and  (14)).  Employing  Equations  (26)^and  {27) 
into_^Eguations  (24)  and  (25),  respectively  together  with  nxnxz=-z, 
and  nxriXT=-T,  one  obtains  the  final  form  of  the  integral  equations  as 

- /HQ2(QlQ‘)Z5(Q)Hj(Q|P)d4=-Hj(Q'  lP)+H^j^(Q'  iP). 

^1 

for  TE^  case  with  Q,Q'  on  Sp  (28) 
and 

^1 

for  TM^  case  with  Q,Q'  on  S-| , (29) 


where  the  following  definitions  have  been  used  in  Equations  (28)  and 
(29). 


'^oz  ^ ^ "oz 


>^6  ^ ^ «ob 


K "b- 


(30a; 30b; 30c) 


Fredholm  integral  equations  of  the  second  kind;  hence,  they  could  be 
solved  via  an  iterative  procedure.  In  particular,  th^  kernels  in  ttie 
above  integral  equations  correspond  to  the  term^  Ho7(q1Q')  and  Hot((51Q') 
in  Equations  (28)  and  (29),  respectively.  Hoz(QiQ'y  represents  the 
z-directed  magnetic  field  at  Q on  the  un_p£rturbed  surface  which  is 
produced  by  a magnetic  line  spurce  z6(jp-pq^l)  at  the  point  Q'  on  the 
same  surface.  Likewise,  Hox(qIQ')  is  the  x-directed  magnetic  fi£ld  at 
0 on  the  unperturbed  surface  due  to  a magnetic  line  dipole  t6(1p-pqi1) 
at  Q'  on  the  same  surface.  Thus,  the  kernel  is  related  to  the  2-D  surface 
Green's  function  for  the  unperturbed  boundary.  The  source  term  Hg[j(Q' iP) 
in  Equation  (28)  represents  the  z-directed  magnetic  field  intensity  induced 
at  Q'  on  th^  jjnperturbed  boundary  S due  to  a unit  amplitude  magnetic  line 
source  zi(lp-pp!)  at  P exterior  to„S.  Similarly,  the  source  term 

Equation  (29)  is  the  x-directed  magnetic  field  intensity 
induced  at  Q'  on  the  unperturbed  boundary  S due  to  a unit  amplitude 
electric  line  source  z5({J-ppi)  at  P exterior  to  S.  Clearly,  in  the 
absence  of  the  perturbation,  i.e.,  when  Zs=0,  the  LHS  of  Equations  (28) 
and  (29)  reduce  identically  to  zero  so  that 
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and 


> 


hJ(q'  Ip)  = jp) 

H^CQ*  iP)  = iP) 


if  Z3  = 0, 


as  would  be  expected.  The  integral  equations  (28)  and  (29)  could 
also  have  been  derived  via  an  application  of  the  compensation  theorem 
for  electromagnetic  fields  [2].  Due  to  the  nature  of  these  equations,  they 
will  be  referred  to  in  this  report  as  perturbation  type  integral  equations. 

The  quantities  Hoz($1Q')  and  h25(Q'1p)  in  Equation  (28),  and 
likewise,  HQx(glQ')  and  H§[j(Q' iP)  in  Equation  (29)  are  assumed  known. 

These  field  quantities  are  known  exactly  for  the  case  of  planar,  and 
circularly  shaped  perfectly-conducting  boundary  S which  are  considered 
in  this  report.  In  the  present  analysis,  the  far  zone  approximations 
to  HCfj(Q'lP)  and  H^[3(Q'jP)  will  be  employed;  furthermore,  since  the 
exact  eigenfunction  series  in  terms  of  the  radial  propagation  repre- 
sentation for  the  fields  pertaining  to  the  perfectly-conducting  circular 
cylinder  geometry  are  poorly  convergent  for  large  cylinders,  the  GTD  repre- 
sentation for  these  fields  is  employed  instead.  Also,  the  GTD  representation 
is  applicable  to  perfectly-conducting  convex  cylinders  of  non-constant  (or 
variable)  curvature  thereby  making  it  more  useful  for  the  present  application. 
The  TE^  and  the  TM^  cases  corresponding  to  the  antenna  problems  of  Figures 
1 and  2 are  next  solved  in  terms  of  their  respective  scattering  problems. 

The  TE^  Problem 


A.  Planar  boundary  case:  The  scattering  problem  from  which  one 
constructs  the  solution  to  the  problem  of  the  radiation  from  an  axial 
slot  in  the  planar  boundary  of  Figure  1 is  analyzed  in  this  section. 

The  geometry  of  this  TEz  scattering  problem,  and  it's  associated 
"unperturbed"  problem  corresponding  to  Figures  4(a)  and  4(b)  are  indi_- 
cated  in  Figures  5(a)  and  5(b),  respectively.  The  source  'r=R=z5 ( !o-pp i ) 
for  the  TE2  scattering  problem  is  shown  in  Figure  5(a),  and  it  is  in  the 
far  zone  of  the  impedance  patch  so  that  the  fields  due  to  R at  Q'  (at  y=0; 
X,x I<L/2J_  ^re  locally  plane.  On  the  other  hand,  the  test  source 
M7=z5( 1p-pq'  j ) is  at  Q'  on  the  unperturbed,  or  perfectly-conducting 
boundary  as  illustrated  in  Figure  5(b).  The  T-direction  indicated  in 
Figures  4(a)  and  4(b)  is  taken  to  be  along  x in  Figures  5(a)  and  5(b), 
respectively  (as  this  corresponds  to  the  increasing  x-coordinate  direction). 
Thus,  Equation  (28)  becomes  in  this  case,  the  following. 


L/2 

= J ‘^^C'5(x-XQ.)-Z^(x)H^^(x!xg,)]H^(x|P);  |x,  and  Xg 
-L/2 


i<L/2 

(32) 
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F (P) 


(a)  PERTURBED  PLANAR 
CONFIGURATION 


(b)  UNPERTURBED  PLANAR 
CONFIGURATION 


Figure  5.  The  scattering  geometry  corresponding  to  the  planar 
perfectly-conducting  surface  with  a surface  imped- 
ance patch,  and  it's  corresponding  unperturbed 
problem. 


in  which  the  coordinates  of  Q'  are  (x=xq';  y=0),  and  the  coordinates 
of  0 are  (x=x;y=0) . 

The  field  Ho2(x1xq')  at  a ooint  Q on  the  unperturbed  (perfectly- 
conducting)  surface  y=0  due  to  % at  Q'  also  on  the  same  surface  is 
easily  shown  to  be 


~r 


(klx-Xg, 


) = 


;x) 


(33) 


(2) 

where  Yq  is  the  free  space  admittance,  ' is  the  cylindrical  Hankel 
function  of  order  zero,  and  k is  the  free  space  wavenumber.  Also, 
^Sb(''<Q‘iP)  '-^hich  is  the  £ield  induced  at  Q'  on  the  unperturbed  surface 
(y=0)  due  tc  the  source  M at  P external  to  S is  given  by 


Q' 


jkxQ,cos4>^ 


(34) 


for  the  case  when  the  source. M is  in  the  far  zone  of  the  imoedance  patch 
(,x|<L/2;  y=0).  The  angle  in  Equation  (34.)  is  illustrated  in  Figure 
5.  Then,  Employing  Equations  (33)  and  (34)  in  (32),  one  may  solve  for 
the  unknown  H'5(x|P)  via  the  conventional  MM  procedure,  provided  Zjix)  is 
assumed  to  be  a known  quantity. 
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In  the  present  MM  solution,  the  unknown  H5(x!P)  is  expanded  in 
terms  of  a sequence  of  pulse  functions. 


kY  — : — N 

hJ(x|P)  Si  - |xl<L/2  and  L-Na, 


The  function  Pq(xq)  is  a unit  amplitude  pulse  function  of  width  A, 
which  is  centered  at  the  point  Xq. 


1 0 i |x-Xj^i> 


Cq  represents  the  unknown  complex  coefficients  in  the  above  expansion 
for  Hg(xlP).  Let  the  point  Q'  (at  y=0;  x=xq')  take  on  the  following 
values. 

Xq.  = Xp  ; p = 1,2,---N,  (37; 

with  the  understanding  that  Xp  corresponds  to  the  mid  point  of  the 
p^l^  pulse  Pp(xp).  This  corresponds  to  point  matching  (or  delta  function 
testing)  at  the  mid  points  of  the  pulse  expansion  functions.  Defining 

jkx„cos^^ 

fp  = e P . (381 


r-  -1 

J dx  z^(x)|^-j2.H(2^kiXp-xl) 


Lno  = - P-V. 

pq  pq 


Xq.A/2 


= I ^ ’P=P 
PP  lo.p^q 


and  incorporating  the  results  of  Equations  (33)  - (39)  into  Equation 
(32)  yields  the  following  matrix  equation  for  Cq. 


[fp]  = [L  ] [C  ] 

P Nxl  PP  MxN  P Nxl 
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It  is  noted  that  the  integration  in  Equation  (39)  is  over  the 
pulse  region,  and  P.V.  denotes  the  principal  value  of  the  integral. 
Finally,  a numerical  inversion  of  the  matrix  [Lpa]  leads  to  a solution 
for  '•He  unknown  Cq. 

B.  Circular  cylinder  case:  The  scattering  problem  from  wnich  one 
constructs  the  solution  to  the  problem  of  the  radiation  from  an  axial 
slot  in  the  circular  cylinder  geometry  of  Figure  2 is  analyzed  in  this 
section.  The  geometry  of  this  TE^  scattering  problem,  and  it's  as- 
sociated "unperturbed"  problem  corresponding  to  Figures  4(a)  and  4(b) 
are  indicated  in  Figures  6(a)  and  6(b)  respectively.  The  source 


(a)  PERTURBED  CYLINDRICAL  (b)  UNPERTURBED  CYLINDRICAL 

CONFIGURATION  CONFIGURATION 


Figure  6.  The  scattering  geometry  corresponding  to  the  circular 
cylinder  oerfectly-conducting  surface  with  a surface 
impedance  patch,  and  it's  corresponding  unperturbed 
problem. 

F’=fT=Z5 ( ;o-7p  1 ) for  this  TE2  scattering  problem  in  Figure  6(a),  and  it 
j_s  in  t^e_far  zone  of  the  cylinder;  whereas,  the  test  source 
Mj=zo ( ip-pqi i ) is  at  Q'  On  the  unperturbed,  or  perfectly-conducting 
cylinder  as  shown  in  Fiqure  6(b).  In  the  present  case,  the  r-direction 
in  Figures  6(a)  and  6(b)  is  chosen  to  corresoond  to  the  p direction 
(as  this  corresponds  to  the  increasing  arc  length  coordinate  a p). 
Hence,  in  this  case.  Equation  (28)  becomes 
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^ob 


(a.$' IP) 


0 


a d<j.[o(i>-V).Z^(«)HQ^(a,<<.la,(?')]Hj(a,^  iP); 

0<<})  ,9 '<24Iq  , (41) 


where  the  coordinates  of  Q'  and  ^ on  the  boundary  S']  are  (a,?')  and 
(a,?),  respectively,  and  9q  is  indicated  in  Figure  6(a). 

The  quantities,  HQ2(a,9 1 a,9 ' ) and  HQ5(a,9’lP)  are  obtained  in  terms 
of  the  GTD  results  which  have  been  developed  elsewhere  [3].  In 
particular,  HQ^(a,9 !a,? ' ) is  given  in  terms  of  H^yp  [3]  as 

CJ 

-j  [a5(t')+jk]dt' 

(a,9|a,9')  ^ C Yq  I Lj(a,9')e  A5(a.9) 

Z— 1 


= Hj^P(a,(fi'  la, 9) 


(42) 


u u u 

where  l2,  A^,  and  aj  are  the  launching,  attachment,  and  attenuation 
coefficients,  respectively  for  Keller's  GTD  surface  ray  mode  fields 
with  modal  index,  4.  The  superscript  "typ”  denotes  the  field  assoiziated 
with  a typical  surface  ray  trajectory  from  Q'  to  Q as  in  Figure^; 


Figure  7.  Ray  paths,  around  the  cylinder  for  calculating 
la,4>' ) or  HQ_(a ,9  ! a ,0  ' ) . 
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one  notes  that  there  are  two  such  ray  paths,  and  HQ2(a,<))  ia.iji' ) is  the 
sum  of  the  fields  associated  with  each  of  these  ray  paths. 

It  is  noted  that  for  numerical  calculations  one  generally  needs 
to  include  only  the  first  couple  of  modes  (corresponding  ^o  ^=1  and 
i=2)  to  obtain  sufficiently  accurate  results  when  Q'  and  (J  are  not  very 
close  to  each  other.  The  higher  order  modes  are  rapidly  attenuated  due 
to  the  fact  that  the  attenuation  constant  increases  significantly 
as  I increases.  For  the  circular  cylinder  case,  is  a constant  for 
a given  radius,  "a",  and  a given  mode  (or  a).  The  quantity  dt'  in 
Equation  (42)  is  the  incremental  arc  length  corresponding  to  the  total 
arc  Q'Q.  The  quantity  H^jj(a,4i ' |P)  is  obtained  in  terms  of  Hoz(pIQ') 
which  is  the  field  at  P due  to  a magnetic  line  source  at  Q'  on  S,  via 
the  relationships  in  Equation  (24)  or  Equation  (26)  and  Equations 
(30a;30b);  furthermore,  HQ2(PiQ')  is  given  in  terms  of  H^|P(PlQ')  [3]  as 


hJ^P(PIQ')  C F 

OZ  ' ' 0 ^ 


jks 


(43) 


where 

Q 

"j  )+Jk]dt' 

F=  I lJ(Q')  e d5(Q); 

i=^ 

if  P is  in  the  shadow  region  of  Q' , 
and 


F = 2 


if  P is  in  the  lit  region  of  Q'. 


(44a) 

(44b) 


u ^ u 

d2(Q)  in  Equation  (44a)  is  the  diffraction  coefficient  for  the 
surface  ray  mode  at  the  point  Q on  the  boundary  S.  The  constant  C 
appearing  in  Equations  (42)  and  (43)  is  given  by  [3] 


C 


k 

/8^ 


(45) 


The  superscript' "typ"  in  Equation  (43)  denotes  the  field  associated  with 
a typical  ray  trajectory  from  Q'  to  P as  in  Figure  8;  one  notes  that  in 
the  lit  region  (see  Figure  8a),  there  are  two  surface  diffracted  ray 
contributions  corresponding  to  ray  paths  Q'QbQaP  and  Q'Q^QbP.  respectively, 
in  addition  to  the  direct  ray  path  Q'P.  The  field  of  the  surface  dif- 
fracted rays  is  given  by  Equation  (43)  with  F as  in  Equation  (44a), 
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(a)  RAYS  FROM  O'  ARRIVING 
AT  P IN  THE  LIT  ZONE 
OF  o' 


RAYS  FROM  Q ARRIVING 
AT  P IN  THE  SHADOW  ZONE 


Figure  8.  Ray  paths  for  calculating  H^^(PlQ')  or  EQ^(PjQ'), 


whereas,  the  direct  ray  field  is  given  by  Equation  (43)  with  F as  in 
Equation  (44b).  The  distance  s in  Equation  (43)  is  si  or  S2  corresponding 
to  QaP  and  QgP,  respectively  for  the  two  surface  diffracted  rays,  and 
s=So=Q'P  for  the  direct  ray.  In  the  shadow  region  (see  Figure  8b),  only 
the  two  surface  diffracted  rays  are  present.  The  fie  Id  Hoz(PiQ')  in  the 

lit  or  shadow  regions  is  thus  the  sum  of  the  fields  Hn|’^(PlQ')  of  all 
the  rays  present  in  Figures  8a  or  8b,  respectively.  The  superscript  h on 
lb,  AO,  and  ot^  refers  to  the  acoustic  hard  type  surface  ray  mode  which 


is  exci 


5, 

ted  in 


a ^ 
the  TE, 


case  under  consideration. 


As  Q'  approaches  (or  vice-versa)  in  Equation  (42),  the  field  point 
lies  in  the  transition  region,  and  the  surface  ray  mode  representation 
in  Equation  (42)  becomes  poorly  convergent;  thus,  it  is  replaced  by  a 
rapidly  convergent  Fock  type  integral  characterizing  the  transition 
region  as  follows: 


JI7 


^2^")  ,-jCt  , 

, .'iv  j\-  e UT 


p:^“(q')  ^ ”2 


where 


f^TT  dt‘  ; ni(t')  = — ^ 


(47a;47b) 


and  Pg(t')  is  the  radius  of  curvature  of  the  convex,  cylindrical 
contour  S at  t‘.  The  quantity  W2(t)  is  the  Fock  type  Airy  function,  and 
W2(t)  is  it's  derivative  [3]. 

Following  an  analysis  which  is  based  on  some  earlier  work  by 
H.  Bremmer,  and  by  Wait  [4],  one  obtains  the  following  approximation  to 
the  Fock  type  integral  of  Equation  (46). 


(a,$la,<^')  'V 


1 


/irj  ,3/2  , .7  ,3 

T-  " ^ ^ 


^ 7/-Trj  9/2  ^ 

+ J . 

with  t = arc  ^Q' . (48) 

Clearly  Equation  (48)  converges  rapidly  for  small  5;  i.e.,  as 
and  it  produces  the  proper  singularity  at  the  source  (at  Q').  Further- 
more, the  first  term  in  the  above  expansion  represents  the  result  for 
the  planar  case  as  given  earlier  in  Equation  (33);  and  the  additional 
terms  in  the  expansion  account  for  the  effects  of  surface  curvature. 

The  transition  region  adjacent  to  the  shadow  boundary  which  is 
associated  with  the  calculation  of  the  field  H^|P(P|Q')  in  Equation  (43) 
also  requires  that  the  GTD  ray  solution  in  terms  of  Equations  (44a) 
and  (44b)  be  modified  within  this  transition  region  as  the  field  point  P 
moves  from  the  shadow  to  the  lit  region  (or  vice  versa)  across  the 
shadow  boundary.  The  simple  ray  descriptions  for  the  fields  fail  within 
transition  regions  across  which  the  different  ray  representations  must 
change  rapidly  but  smoothly.  A uniform  result  for  the  transition  region 
which  reduces  to  the  proper  ray  field  descriptions  in  the  lit  and 
shadow  zones  outside  the  transition  region  (as  required  by  Equations  (44a) 
and  (44b))  is  presented  in  [3];  however,  a more  convenient  representation 
based  on  a direct  extension  of  Ivanov's  work  [5]  is  presented  below. 
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jkdt' 


F = g(5)e 


Q* 


wi  th 


py^(Q') 


if  P is  in  the  shadowed  part  of  the 
transition  region  (i.e.,  for  5>0), 

(49a) 


5 = 


m(t' ) 


dt' 


Likewise, 


F = g(0 


with 


5 = - m(Q' ) |n-s 


(49b) 


if  P is  in  the  lit  part  of  the 
transition  region  (i.e.,  for  C<0), 


(50a) 

(50b) 


The  unit  vectors  n and  s are  indicated  in  Figure  8a.  The  acoustic  hard 
Fock  function,  g(5)  is  tabulated,  and  it  is  defined  in  [3].  It  may  be 
verified  that  Equation  (49a)  reduces  uniformly  to  Equation  (44a)  in  the 
deep  shadow,  whereas  Equation  (50a)  reduces  uniformly  to  Equation  (44b) 
in  the  deep  lit  region.  Employing  the  origin  as  the  phase  reference 
for  the  circular  cylinder  case  under  consideration.  Equation  (43)  may 
be  rewritten  in  the  far  zone  as 


hH^(pIQ')  C F 
OZ  ' ' 0 f- 

v D 


-jkp 


(51a) 


where  the  phase  factor  is  easily  seen  to  be 

^ ka  cos[i{i ' -({i’ ] ; if  P is  in  the  lit  region  of  Q* 


0 


; if  P is  in  the  shadow  region  of  Q' 


(51b) 


The  angles  and  are  indicated  in  Figures  8a  and  8b 
m Hk(a, 

pulse  functions;  thus. 


The  unknown  Hk(a,(flP)  is  expanded  as  before  in  terms  of  a set  of 
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(52) 


r 


-jkp  N 

HK(CfiP)  ^ C Y„  • I C P J ; 


0 J~ 

/p  q=l 


q q'^q' 


where  the  unit  amplitude  pulse  function  Pq('<>q) 
centered  at  the  point  (a.ifg). 


1 i 
0 ; 


0<|«  |<2$^ 
and  2(tiQ=NA 

of  angular  width  a is 


(53) 


In  order  to  evaluate  the  unknown  coefficients  Cq  in  the  expansion  of  : 

Equation  (52),  one  employs  point  matching  (or  delta  function  testing)  i 

at  the  center  of  each  pulse;  i.e.,  one  lets  the  point  Q'  at  (a.^i')  ' 

take  on  the  following  values,  i 

I 

I 

= <(>p  ; p = 1 ,2,*  • - N , (54)  i 

with  the  understanding  that  (a,i>p)  is  the  center  of  the  p^^  pulse 
Pp('^p)  • 

As  done  before  for  the  planar  case,  one  may  incorporate  the 
results  of  Equations  (42)  - (54)  into  Equation  (41)  to  arrive  at  a 
matrix  equation  analogous  to  that  in  Equation  (40)  from  which  the  unknown 
Cq  may  be  found  via  a numerical  matrix  inversion  procedure.  It  is  noted 
that  the  arc  length  integrals  in  Equations  (42),  (44a),  (46),  (47a), 

(49a)  and  (49b)  are  trivially  integrated  in  closed  form  for  the  circular 
cylinder  case;  however,  the  general  forms  of  these  integrals  have  been 
presented  to  indicate  that  the  above  method  of  solution  is  also  applicable 
to  convex  cylinders  of  nonconstant  curvature. 

The  TM^  Problem 


A.  Planar  boundary  case:  The  scattering  problem  from  which  one 
constructs  the  solution  to  the  problem  of  the  radiation  from  a circum- 
ferential slot  in  the  planar  boundary  of  Figure  1 is  analyzed  in  this 
section.  The  geometry  of  this  TMz  scattering  problem,  and  it's  associated 
"unperturbed"  problem  corresponding  to  Figures  4^  and  4b  are  indicated 
j_n  Figur£S_5a  and  5b,  respectively.  The  source  F is  now  of  the  type 
F_=T=zq(  |p2P£|  ) in  the  TM2  case;  whereas,  the  test  source  ITj  is  given  by 
Mj  = xq(Ip-pqiI)  at  Q*  on  the  unperturbed^ (or  perfectly-conducting) 
surface  with  the  t direction  being  along  x in  this  case.  The  corresponding 
integral  equation  for  the  TM^  case  becomes  via  Equation  (29): 
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^ob^^Q' 


L/2 

-1/2 


dx[6(x-XQ.)-Z^(x)H^^(x!xQ.)]H®(x!P); 

|xi<L/2 

1xq.1<L/2  , 


(55) 


in  which  the  coordinates  of  Q'  are  (x=xq'  ; y=0),  and  the  coordinates 
of  Q are  (x=x;y=0). 

The  field  HoT(x|xq')  at  on  the  unperturbed  surface  due  to  Mj 
at  Q'  also  on  the  same  surface  is  easily  shown  to  be 


kY 


,(2) 


0 


H„,(xlXjj,)  = - 'kix-x  ° N„.(Xni  lx) 


OT 


ot'''Q' 


(56a) 


in  which  is  the  cylindrical  Hankel  function  of  order  one.  It  is 
noted  that  Hox(xlxq')  is  more  singular  than  HQ2(xixq')  of  Equation  (33) 
for  the  TE^  case;  consequently,  an  alternative  representation  for 
^^OT(xixq')  which  will  be  employed  for  xq'-*-x  in  the  present  analysis  is 


H,^(Xq.(x)=.  ° 


kY 

( 

T" 


HQ^^(k|x-XQ,  I)  + ^ HjJ^^(k|x-Xg,  I) 

_ _ 


(56b) 


The  field  H®. (xqi  iP)  due  to  J which  is  induced  at  Q'  on  the 
unperturbed  surface  is  given  by 


sin,'  e 


jkXg,cos4) 


(57) 


when  the  source  J is  in  the  far  zone  of  tht  surface  impedance  patch 
region.  The  angle,  is  illustrated  in  Figure  5.  One  may  incorporate 
Equations  (56b)  and  (57)  into  Equation  (55)  and  solve  for  H^(xlP)  via 
the  MM  procedure;  thus,  one  may  expand  H^(xlP)  in  terms  of  a sequence 
of  pulse  functions  as 


i '^q  |xhl-/2,  and  L=Na.  (58) 


1 irko 


,-jkp 


N 


q=l 


q q'  q' 


The  function  Pq(xq)  is  defined  previously  in  Equation  (36);  and  Cq  are 
the  unknown  coefficients  which  are  to  be  determined. 
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Let 


Also, 


X *i/2 

P.V.  J dx  Zj(x)  H<^'(kix-Xp!) 
X -A/2 

q 


(59) 


3x5. 


9q(^Q') 


^ 3n(Xn) 

-3 — E_- 


2ga(x  ) + g (x") 

9-£—  ; at  Xq.=Xp  (60) 


XQ,=Xp 


where 


h s jIA/2  ; wi  th  , 


(61) 


and 


x"  = Xp  ± h (specifically,  Xp  is  defined  in  (63)  below).  (62) 


A finite  difference  approximation  has  been  employed  in  Equation  (60),  and 
the  finite  difference  width,  2h  depends  on  the  parameter  l as  indicated 
in  Equation  (61).  The  quantity  a in  Equation  (61)  of  course  corresponds 
to  the  width  of  the  expansion  pulse  Pq(xq). 

In  the  present  case,  the  choice  of  i is  critical,  and  the  value  of 
^ chosen  is  approximately  23/60;  more  will  be  indicated  on  this  choice 
of  i for  obtaining  stable  numerical  results  in  Section  III.  The  unknown 
Cq  may  be  found  by  solving  a set  of  simultaneous  linear  equations 
obtained  by  allowing  the  point  Q*  (at  x=xqi,  y=0)  to  take  on  the  values 


Xq i”Xp  , p—  1,2,*""N, 


(63) 


which  corresponds  to  point  matching  at  the  center  of  the  pulses.  Defining 


. jkx  coS(>^ 
ap  = siniji  e ^ 


and 


kY 

N =6  + 

pq  pq  ~ 

Xoi  =X 

Q P 
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T 


(64) 


(65) 


one  may  express  the  set  of  equations  for  Cq  compactly  in  a matrix 
notation  in  terms  of  ap  and  Npq  as 


[N  ] 
pq 


NxN 


(66) 


Thus,  Cq  may  be  obtained  via  a numerical  matrix  inversion  of  [Npq], 

B.  Circular  cylinder  case:  The  scattering  problem  from  which  one 
constructs  the  solution  to  the'  problem  of  the  radiation  from  a circum- 
ferential slot  in  the  circular  cylinder  geometry  of  Figure  2 is 
analyzed  in  this  section.  The  geometry  of  this  TM2  scattering  problem, 
and  it's  associated  "unperturbed"  problem  corresponding  to  Figures  4a  and 
4b  are  indijcated  in  Figures  6a  and  6b  respectively.  The  source 
r=T=ziS(  |p-pp  I ) for  this  TMj  scattering  problem  in  Figure  6a,  and  it_i£ 
in  the  far  zone  of  the  cylinder;  whereas,  the  test  source  % = i'5(!p-PQ'|) 
is  at  Q'  on  the  unperturbed  (or  perfectly-conducting)  cylinder  as  shown 
in  Figure  6b.  The  t-direction  has  of  course  been  chosen  to  correspond 
to  the  ^-direction.  Thus,  the  TM2  integral  equation  for  this  problem 
becomes  via  Equation  (29)  the  following. 


Hob(a«'^'lP)  =j  ad45(^-V)-Zj(^)HQ^(a,<p;a,«')]H®(a,<^lP); 

0 

0<<j)  ,ii '<241^,  (67) 

where  the  coordinates  of  Q'  and  0 on  the  boundary  S,  are  (a,$')  and 
(a,())),  respectively  as  in  the  TE2  case  (see  Equation  (41));  and  pq  is 
as  in  Figure  6a. 


Following  the  procedure  for  solving  the  TE2  cylinder,  one  begins 
by  incorporating  the  GTD  approximations  for  Hox(a,i>  la,i) ' ) and 
H®jj(a,(|)' j P) . Thus,  H^'yP(a,(ti  ja,p' ) where  the  superscript  "typ"  has  the 
same  meaning  as  in  Equation  (42)  for  the  TE2  case  is  given  by  [3] 


HjyP(a,$la,i.')=CYQ  lL^(a,<fr')e 

= HjyP(a,<.'la,<^) 


[ct^(t'  )+jk]dt' 


(68) 


c s s 

Aj^  and  are  now  the  launching,  attachment,  and  attenuation 
coefficients  for  the  4^-^^  GTD  surface  ray  mode.  Again,  an  inclusion  of 
only  the  first  couple  of  modes  ;is  sufficient  for  obtaining  accurate 
results  when  the  points  Q'  and  Q are  not  too  close  to  each  other. 
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HQ^(a,(i)la,i)>' ) is  the  sum  of  the  fields  H^'^Passociated  with  each  of  the 
two  ray  paths  illustrated  in  Figure  7. 

The  quantity  H®|j(a,(ji ' 1 P)  is  obtained  in  terms  of  EQ^(PiQ')  which 
is  the  field  at  P due  to  a magnetic  line  dipole  source  at  Q'  on  S via 
the  relationships  in  Equations  (25)  or  (27)  and  (30a;30b);  in  particular 
EiyP(PlQ')  is  given  by  [3]: 


eJ^P(pIQ')  'V  c G 


(69) 


where  the  meaning  of  the  superscript  "typ"  and  the  distance  s is  the 
same  as  in  Equation  (43),  and  where 

q. 

-J  Ca^(t' )+jk]dt' 

G = I L^(q')  e ^ D^(Q);  if  P is  in  the 

^=1  shadow  region  of  Q',  (70a) 


and 


G = 2 n*s  ; if  P is  in  the  lit  region  of  Q'. 


(70b) 


The  d|(Q)  in  Equation  (70a)  is  the  diffraction  coefficient  of  the 
surface  ray  mode  at  Q.  The  constant  C in  Equation  (69)  is  given 
previously  in  Equation  (45),  and  the  superscript  s on  L|,  a|,  D|  and 
a|  refers  to  the  acoustic  soft  type  surface  ray  mode  which  is  excited 
in  the  TM,  case  under  consideration.  The  ray  trajectories  associated 
with  E^P(PIQ')  in  Equation  (69)  are  indicated  in  Figure  8,  and  Eq^  is 
the  sum  of  the  fields  Eq^P  which  are  associated  with  the  rays  shown  in 
this  figure. 


As  in  the  TEz  case,  when  Q'  and  are  sufficiently  close,  one  is 
within  the  transition  region,  and  Equation  (68)  must  be  replaced  by  a 
more  rapidly  converging  Fock  type  integral  as 


HQ^P(a,<(.|a,$')^-CYQe 


w^(t) 

wjrry 


(71) 


where  the  quantities,  5,  Pg,  m,  W2 
as  in  Equation  (46).  An  extension 
yields  the  following  approximation 


and  w^  appearing  above  are  the  same 
of  the  analysis  in  [4]  to  the  TMz  case 
to  Equation  (71 ) as 
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(72a) 


krt 


iy|I^3/2^3,5^^9/2^... 


The  above  approximation  is  analogous  to  that  in  Equation  (48)  for 
the  TEz  case;  and  as  in  Equation  (48)  it  is  seen  that  the  first  term  of 
the  expansion  in  Equation  (72)  represents  the  flat  or  planar  surface 
result  corresponding  to  Equation  (56a),  with  the  additional  terms 
representing  the  correction  to  the  flat  surface  result  arising  from  the 
presence  of  finite  curvature  of  the  cylinder.  Following  the  alternative 
representation  for  Equation  (56a)  given  in  Equation  (56b),  one  may  also 
rewrite  Equation  (72a)  for  the  circular  cylinder  as: 


kY 

HQ'^*^(a,(j) ' |a,(^)''^ — ^ 


H^^^klti)  + ^ L-^H^^^klti) 

° (ka)'^  ° 


1 {H  r3/2  ^ j5  ^3  ^ 5 r-T  ^9/2  * 


(72b) 


where  lt|  = a|(^-<()'l  in  the  circular  cylinder  case.  The  above  could  be 
generalized  to  arbitrary  convex  cylinders  using  the  local  properties 
of  asymptotic  high  frequency  results.  The  result  in  Equation  (72b) 
is  useful  for  large  ka  and  small  lt|  (i.e.,  for  small  c). 

One  must  also  modify  the  function  G in  Equation  (69)  as  the  point  P 
moves  within  the  transition  region  across  the  shadow  boundary  as  done 
for  the  function  F in  Equations  (49)  and  (50)  for  the  TEz  case.  Based  on 
the  work  in  [5],  one  obtains 


G 


jkdt' 


(Q) 

(O')  ’ 


with  5 as  in  Equation  (49b);  and. 


if  P is  in  the  shadowed 
part  of  the  transition 
region  (i.e.,  for  5=*0) , (73) 


G = g(5)  e'3  / ; if  P is  in  the  lit  part 

' of  the  transition  region 

(i.e.,  for  5<0) , (74) 

with  5 as  ^n  Equation  (50b),  respectively.  The  acoustic  soft  Fock 
function,  g(c)  is  tabulated,  and  it  is  defined  in  [3]. 
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It  may  be  verified  that  Equation  (73)  reduces  to  Equation  (70a)  in 
the  deep  shadow,  whereas  Equation  (74)  reduces  to  Equation  (70b)  in  the 
deep  lit  region  thereby  providing  a uniform  approximation  to  G within 
the  transition  region  across  the  shadow  boundary.  Using  the  origin 

p-Jks  , ^ .-jkP 

the  phase  reference,  one  replaces  in  Equation  (69)  by 


as 


where  is  defined  in  Equation 


(51b). 


1 

i 


The  unknown  H^(a,itlP)  is  then  expanded  in  terms  of  a set  of  pulse 
functions  as 


,-jkp 


h®(q!p)  4 c I c p (♦  ) 

✓p  q=i  -I  '1  M 


(75) 


analogous  to  that  in  Equation  (52)  for  the  TE^  case.  For  Itl  or  C'Kl, 
the  approximation  in  Equation  (72b)  is  employed  instead  of  Equation  (72a), 

and  the  5-  derivative  is  approximated  via  the  finite  difference 

3(i  ' 

scheme  similar  to  that  indicated  in  Equation  (60)  for  the  planar  TM^ 
case. 

Then  employing  point  matching  (or  impulse  testing)  at  4i'='jip, 
p=l,2,"*-N  as  in  the  TE2  case,  one  arrives  at  a matrix  equation  of  the 
form 


P Nxl 


(76) 


analogous  to  the  planar  TM^  case  from  which  the  unknown  Cq  may  be  found 
by  a numerical  inversion  of  [Npq].  The  explicit  forms  or  a^  and  Npq 
are  not  presented  for  the  sake  of  brevity,  but  they  are  easily  found 
from  the  results  for  Hqt  and  Eqt  presente(j  above.  The  finite  difference 
width  is  chosen  to  be  the  same  as  that  determined  for  the  planar  TMz 
case. 
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III.  NUMERICAL  RESULTS  AND  DISCUSSION 


Having  obtained  the  solutions  to  the  surface  field  induced  in  the 
problem  of  scattering  of  electromagnetic  waves  by  planar  and  cylindrical 
perfectly-conducting  surfaces  which  are  covered  by  a surface  impedance 
patch  of  finite  extent,  one  may  then  directly  calculate  the  radiation 
pattern  of  a thin  slot  placed  within  the  surface  impedance  region  of 
these  configurations  via  the  reciprocity  principle  as  indicated  in  the 
previous  section. 

Basically,  the  radiation  pattern  calculation  amounts  to  finding  the 
amplitude  Cg  of  the  expansion  pulse  representing  the  value  of  the 
induced  surface  magnetic  field  in  the  scattering  problems  such  that  this 
qth  pulse  is  centered  at  the  slot  location  in  the  corresponding  antenna 
problem.  The  radiation  pattern  of  this  slot  in  the  if''  direction 
(see  Figures  5 and_6)  is  th^n  obtained  by  calculating  the  value  of  Cg 
due  to  the  source  F at  (p.if^)  in  the  scattering  problem.  In  the  present 
calculations,  Cg  is  calculated  for  every  two  degree  increment  in 

The  moment  method  solution  of  the  perturbation  integral  equations 
is  carried  out  using  pulse  expansion  functions,  and  impulse  (or  delta) 
testing  functions;  this  choice  of  expansion  and  test  functions  has  been 
primarily  selected  for  simplicity.  As  mentioned  in  the  previous  section, 
the  TM^  solutions  presented  some  numerical  difficulties  due  to  the  more 
singular  nature  of  the  kernel. 

While  pulse  expansion  functions  and  piecewise  sinusoidal  test 
functions  were  initially  employed  to  circumvent  this  difficulty,  it 
turned  out  that  the  numerical  results  were  sensitive  to  the  width  of 
the  test  function  in  this  case.  On  the  other  hand,  the  finite  difference 
approximation  employed  together  with  pulse  expansion  and  impulse  test 
functions  led  to  results  which  were  sensitive  to  the  choice  of  the  finite 
difference  width.  Results  sensitive  to  the  width  of  test  functions, 
and  also  to  finite  difference  intervals  have  been  encountered  elsewhere 
in  the  MM  solutions  to  wire  antenna  problems  [6,7]. 

The  advantages  gained  by  employing  the  perturbation  type  integral 
equation  which  requires  the  unknown  to  be  solved  only  over  the  surface 
impedance  patch  may  thus  appear  to  be  negated  for  the  TMz  case  due  to 
the  occurrence  of  a more  singular  kernel  than  that  which  is  present  in 
the  conventional  electric  field  integral  equation  (EFIE). 

The  EFIE  is  a special  case  of  the  reaction  integral  equation  [8]. 
However,  it  is  shown  that  once  the  proper  finite  difference  width  is 
chosen  (this  width  is  then  the  same  for  all  the  TMz  cases  considered 
herein),  the  results  for  the  TMz  case  are  indeed  stable  and  can  be 
simply  computed.  In  the  present  analysis,  the  finite  difference  width 
was  varied  until  it  yielded  a value  for  a certain  parameter  in  the 
problem  which  agreed  closely  with  that  obtained  via  an  accurate, 
independent  procedure.  One  recalls  that  the  finite  difference  scheme 
is  required  in  the  TMz  case  to  approximate  the  second  derivative  in  the 
kernel  of  Equation  (56b)  or  Equation  (72b). 
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In  actual  numerical  computations,  the  representations  for  the 
kernels  in  Equations  (56b)  and  (72b)  corresponding  to  the  planar  and 
cylindrical  cases  were  employed  only  for  computing  the  elements  Nn, 

Ni2  and  N-i3  (note:  Ni]  = N22  = 1^33  = •••  = Npp)  in  Equations  (66)  and 
(76),  respectively,  whereas,  the  representations  in  Equations  (56a)  and 
(72a)  were  employed  for  calculating  all  other  Npq  elements  with 
q^+3.  The  Npq  for  q<p  are  then  obtained  from  Npq  for  q>p  via 
symmetry  since  Npq  = Nqp. 

The  results  of  the  radiation  pattern  calculation  of  slots  in  the 
planar  and  cylindrical  configurations  of  Figures  1 and  2 are  indicated 
in  Figures  9-22  for  the  axial  slot  (or  TEz)  case,  whereas  those  for  the 
circumferential  slot  (or  TMz)  case  are  indicated  in  Figures  23-29,  I 

respectively.  j 

i 

These  MM-GTD  results  are  compared  with  numerical  results  obtained  i 

by  other  independent  methods  which  are  based  on  either  the  complete  GTD  j 

approach  [11]  or  the  complete  MM  solution  of  the  conventional  reaction  j 

integral  equations  for  the  same  problems.  Also,  an  exact  eigenfunction  1 

series  expansion  is  available  for  calculating  the  radiation  from  slots  j 

in  circular  cylinders  which  are  completely  covered  with  an  impedance  ' 

surface.  Such  an  expansion  has  been  computed  for  a few  TEz  cases  to  ] 

provide  a comparison  with  the  corresponding  MM-GTD  results  as  indicated  j 

in  Figures  12  and  13.  Indeed,  the  agreement  between  the  results  based  | 

on  the  present  W-GTD  method  and  by  other  independent  methods  is  in 
general  quite  good.  Unless  otherwise  specified,  the  patterns  are  normalized  ] 

to  0 dB  at  their  peak  values.  Figure  17  indicates  the  effect  of  changing  | 

the  value  of  the  surface  impedance  patch  on  a cylinder  from  an  inductive  j 

to  a capacitive  type  reactance  for  the  axial  slot  (or  the  TE^)  case. 

It  is  noted  that  a capacitive  reactance  will  not  support  an  Elliot  type 
mode  in  the  TEz  problem,  and  therefore,  the  gain  near  the  horizon  would 
be  reduced  in  this  case. 

Radiation  pattern  calculations  via  the  MM-GTD  method  are  included 
to  illustrate  the  effects  of  changing  the  size  of  the  surface  impedance 
patch,  changing  the  value  of  the  surface  impedance,  and  also  changing 
the  size  of  the  circular  cylinder.  In  general,  it  is  seen  that  in- 
creasing the  size  of  the  impedance  patch  increases  the  number  of  pattern 
ripples  in  the  lit  region,  i.e.,  in  the  region  of  space  where  the  slot  is 
directly  visible-,  whereas, increasing  the  size  of  the  cylinder  tends  to 
decrease  the  level  of  the  radiation  pattern  in  the  deep  shadow  region 
of  the  slot.  Over  the  range  of  impedance  values  considered  in  the 
present  calculations  one  notes  that  increasing  the  magnitude  of  the 
surface  impedance  increases  the  size  of  pattern  ripple  in  the  lit  region 
for  the  TEz  case,  whereas  it  does  the  opposite  for  the  TM2  case.  A means 
to  reduce  the  level  of  the  pattern  ripple  in  the  lit  region  is  to 
appropriately  taper  the  surface  impedance  as  indicated  in  Figure  16. 


31 


s? 

In  conclusion,  a perturbation  type  integral  equation  has  been 
employed  for  solving  the  antenna  problems  illustrated  in  Figures  1 and 
2 via  a combination  of  fli  and  GTD  approaches.  This  procedure  requires 
the  unknown  to  be  evaluated  via  the  MM  procedure  only  over  the  surface 
impedance  patch  region,  hence  it  is  different  from  the  MM  solution  of 
conventional  integral  equations  (for  the  same  problem)  which  employ  the 
free  space  Green's  function,  as  the  latter  requires  the  unknown  to  be 
evaluated  not  only  over  the  surface  impedance  patch  region,  but  also 
over  the  rest  of  the  perfectly  conducting  region.  The  systematic  use  of 
the  GTD  provides  simplifications  in  the  MM  calculations,  and  it  allows 
one  to  treat  convex  cylinders  with  variable  curvature  without  much 
difficulty.  Thus,  the  present  approach  can  treat  electrically  large 
perfectly-conducting  cylinders  (with  impedance  patches)  far  more  ef- 
ficiently than  by  the  use  of  the  conventional  integral  equations.  Further- 
more, the  GTD  approximations  incorporated  in  the  kernel  and  source  terms 
of  the  perturbation  integral  equations  become  increasingly  accurate  as 
the  frequency  increases,  and  in  general  these  GTD  approximations  are 
quite  accurate  even  for  ka  as  small  as  6. 

Further  work  needs  to  be  done  to  arrive  at  a more  satisfactory 
approach  for  obtaining  a proper  choice  of  the  finite  difference  width 
in  the  circumferential  slot  (or  TM^)  case.  For  example,  one  might 
investigate  how  well  the  appropriate  boundary  conditions  are  met  on  the 
perfectly-conducting  portion  as  a function  of  the  finite  difference 
width,  or  investigate  some  other  simpler  criteria  determining  this 
finite  difference  width. 

With  the  present  solutions  one  could  now  begin  to  investigate  the  CP 
characteristics  of  crossed  slots  in  the  surface  impedance  loaded  circular 
cylinders  as  a function  of  the  different  parameters  involved.  In 
addition,  it  would  be  of  interest  to  include  loss  in  the  surface  impedance 
to  study  how  that  would  affect  the  pattern  shape  and  CP  characteristics 
of  crossed  slots  in  cylinders.  The  extension  of  this  analysis  to  include 
airplane  wing  effects  in  the  roll  plane  can  be  carried  out  via  GTD  in 
which  the  fields  incident  on  the  wings  may  be  calculated  via  the  present 
MM-GTD  approach,  while  the  reflections  from  the  wing  surface  and  the 
diffractions  from  the  wing  tips  due  to  the  fields  incident  on  the  wings 
may  be  easily  and  simply  calculated  via  the  ray  reflection  and  diffraction 
coefficients  as  in  [10]. 
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dtliation  pattern  of  a magnetic  line  source  on  a perfectly-conducting  planar 
urfdce  wliicli  is  covered  wi  tit  an  impedance  surface  patch. 
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Figure  13.  Radiation  pattern  of  a magnetic  line  source  on  a 
perfectly-conducting  circular  cylinder  whicn  is 
covered  with  an  impedance  surface  oaten. 


37 


AXIAL  SLOT 
OR  TE2  CAS 

270* 


Figure  14.  Radiation  pattern  of  a 'Magnetic  line  source  on  a 
perfectly-conducting  circular  cylinder  which  is 
covered  with  an  inpedance  surface  patch. 


Figure  21.  Radiation  pattern  of  a magnetic  line  source  on  a 
perfectly-conducting  circular  cylinder  which  is 
covered  with  an  impedance  surface  patch. 
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Figure  22.  Radiation  pattern  of  a magnetic  line  source  on  a 
perfectly-conducting  circular  cylinder  which  is 
covered  with  an  impedance  surface  patch. 
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igure  26.  Radiation  pattern  of  a magnetic  line  dioole  source 
on  a perfectly-conducting  circular  cylinder  which 
is  covered  with  an  impedance  surface  patch. 
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figure  27.  Radiation  pattern  of  a magnetic  line  dipole  source 
on  a perfectly-conducting  circular  cylinder  which 
is  covered  with  an  impedance  surface  patcli. 
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Figure  2H.  Kdclidtiun  pattern  of  a magnetic  line  dipole  source 
on  d perfectly-conducting  circular  cylinder  which 
is  covered  witlt  an  impedance  surface  patch. 
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METRIC  SYSTEM 


1 


BASE  UNITS: 

Quantity 


Unit. 


length 

mass 

time 

electric  current 
thermodynamic  temperature 
amount  of  substance 
luminous  intensity 

SUPPLEMENTARY  UNITS: 

plane  angle 
solid  angle 

DERIVED  UNITS: 
Acceleration 

activity  (of  a radioactive  source) 

angular  acceleration 

angular  velocity 

area 

density 

electric  capacitance 

electrical  conductance 

electric  field  strength 

electric  inductance 

electric  potential  difference 

electric  resistance 

electromotive  force 

energy 

entropy 

force 

frequency 

illuminance 

luminance 

luminous  flux 

magnetic  field  strength 

magnetic  flux 

magnetic  flux  density 

magnetomotive  force 

power 

pressure 

quantity  of  electricity 
quantity  of  heat 
radiant  intensity 
specific  heat 
stress 

thermal  conductivity 
velfK.ity 

viscosity,  dynamic 

viscosity,  kinematic 

voltage 

volume 

wavenumber 

work 


metre 

kilogram 

second 

ampere 

kelvin 

mole 

candela 


radian 

steradian 


metre  per  second  squared 

disintegration  per  second 

radian  per  second  squared 

radian  per  second 

square  metre 

kilogram  per  cubic  metre 

farad 

siemens 

volt  per  metre 

henry 

volt 

ohm 

volt 

joule 

joule  per  kelvin 

newton 

hertz 

lux 

candela  per  square  metre 
lumen 

ampere  per  metre 

weber 

tesla 

ampere 

watt 

pascal 

coulomb 

joule 

watt  per  steradian 

joule  per  kilogram-kelvin 

pascal 

watt  per  metre-kelvin 
metre  per  second 
pascal-second 
square  metre  per  second 
volt 

cubic  metre 
reciprocal  metre 
joule 


SI  PREFIXES: 

Multiplication  Factors 


SI_Symbol 

Formula  _ 

m 

>ig 

s 

A 

K 

mol 

cd 

rad 

sr 

m/s 

(disintegration  )/s 

rad/s 

rad/s 

m 

kg/m 

F 

A-sA/ 

S 

m 

V/m 

H 

V-s/A 

V 

W/A 

V/A 

V 

W/A 

1 

N-m 

]IK 

N 

kg-nvs 

Hz 

(cycle  )/s 

lx 

Im/m 

cd/m 

Im 

cd-sr 

A/m 

Wb 

V-s 

T 

Wb/m 

A 

W 

|/s 

Pa 

N/m 

C 

A-s 

1 

N-m 

Wsr 

I/kg-K 

Pa 

N/m 

W/m-K 

m/s 

Pa-s 

m/s 

V 

W/A 

m 

(wave)/m 

1 

N-m 

Prefix 

SI  Symh«>l 

1 OOO  000  000  000 

:= 

10” 

lera 

T 

1 000  000  000 

s 

10’ 

K'ga 

c; 

M 

1 000  000 

s 

10* 

mega 

1 000 

a 

\0‘ 

kilo 

k 

100 

a 

10' 

hecto* 

h 

10 

10' 

dnka* 

da 

0.1 

10-' 

dec!’ 

d 

0 01 

a 

10-' 

cent)* 

r 

0 001 

10-* 

milll 

m 

0.000  001 

10“‘ 

micro 

M 

0 000  000  001 

10-’ 

nano 

n 

0.000  000  000  001 

10-'' 

pico 

P 

0 000  000  000  000  001 

10-'* 

fomlo 

f 

0 000  000  000  non  ooo  ooi 

in-" 

alto 

a 

To  be  avoided  where  possibli 


